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Abstract 

We show that measurements of finite duration performed on an open two-state system can pro- 
tect the initial state from a phase-noisy environment, provided the measured observable does not 
commute with the perturbing interaction. When the measured observable commutes with the envi- 
ronmental interaction, the finite-duration measurement accelerates the rate of decoherence induced 
by the phase noise. For the description of the measurement of an observable that is incompatible 
with the interaction between system and environment, we have found an approximate analytical 
expression, valid at zero temperature and weak coupling with the measuring device. We have 
tested the validity of the analytical predictions against an exact numerical approach, based on 
the superoperator-splitting method, that confirms the protection of the initial state of the system. 
When the coupling between the system and the measuring apparatus increases beyond the range 
of validity of the analytical approximation, the initial state is still protected by the finite-time 
measurement, according with the exact numerical calculations. 
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I. INTRODUCTION 



The problem of measurement is fundamental to quantum theory lH5|. According to 
;he standard interpretation , established after the pioneer contributions by Born 

3), Dirac and von Neumann j3j|, the postulated collapse of the wave function 
implies an instantaneous measurement process. However, the de facto meaning of a quantum 
measurement is the result of a physical interaction between the system being measured and 
;he measuring apparatus. Such an interaction can be described by a suitable formalism 
9), [loj|, which, under the hypothesis that the measurement process is irreversible, leads to 



the Lindblad equation 



11 



13] . This equation describes the evolution of a measurement 



during a finite time interval by the stochastic term - Lindbladian. At last, the information 
of the probabilities associated to every possible result is acquired by the conventional way 



14 



3- 



16 



,ll2||, we have been 



with the density operator diagonal elements - populations 

Based on superoperator algebra and Nakajima-Zwanzig projectors 
able to describe the dynamics of an arbitrary measurement that occurs during a finite time 
interval, while the system being measured interacts with the rest of the universe and, due to 



the consequent environmentally- induced noise, undergoes decoherence [18|. The assumption 
that the interaction of the measuring apparatus with the system is Markovian justifies a 
Lindbladian approach. However, to treat the noise introduced by the fact that, during 
the finite-duration measurement, the system is perturbed by the environment, a Markovian 
approximation is too restrictive, since non-Markovian noise effects can become non-negligible 



for time scales in which coarse graining is not a good approximation |19H24|. Instead, to 
develop a formalism able to include non-Markovian effects, we have used a Redfield approach 
to the dynamical description of the interaction between the system and its environment, 
excluding the apparatus. After tracing out the degrees of freedom introduced to describe the 
non-Markovian noise, the resulting master equation in the Born approximation (not a Born- 
Markov approximation), referred only to the system, can be used to investigate the effects of 
the environmental noise on the measurement dynamics. We do not use additional Lindblad 
channels to describe the environmental noise because such a formulation would preclude any 
attempt to treat non-Markovian effects. Since we are interested in describing the effects of 
the noise during the finite-time measurement, turning off the environmental perturbation 
during the action of the measuring apparatus is not a valid approximation. Therefore, our 
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hybrid approach is the most economical dynamical description of a measurement that is 
simultaneous to non-Markovian environmental perturbations. 

The Lindblad equation can be used to describe Markovian and non-Markovian environ 



12 



mental effects, with applications ranging from decoherence and dissipation analyses 
to measurement processes particularly when we consider that the measurements have 



finite durations 



151 ] . It follows as a mathematical consequence of a semi-group dynam- 



ics to describe irreversibility in quantum mechanics ll|. Therefore, it can only be used 
to describe Markovian processes, unless the number of dynamical variables describing the 



system be enlarged 



20], causing computational overhead that we must avoid. Incidentally, 
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we point out that the extension of the Lindblad equation to non-Markovian systems 
can be applied to non-semigroup dynamics as well, with time-dependent Hamiltonian and 
Lindbladian operators. 

In the present paper, we analyze finite-time measurements that commute or do not com- 
mute with the interaction Hamiltonian. We study a two-state system interacting with a 
bath of harmonic oscillators, that emulates a phase- noisy environment. The approximate 
analytical solution agrees with the numerical results of the superoperator-splitting method 



25| for weak coupling between the system and its environment. We find that the finite- 



time measurement can protect the measured state, if the observable being measured does 
not commute with the Hamiltonian describing the noisy environmental interaction with the 
system but, when the measurement commutes with the interaction, the environment only 
increases the coherence decay. For strong coupling between the system and the environment, 
our analytical solution fails, as shown by a thorough comparison of its predictions against 
the corresponding exact numerical results. In the case of strong coupling, the effect pre- 
dicted by the numerical calculations is a more intense error as the coupling with environment 
increases. 

Subsequent researchers may use the method presented here to study the effects of varying 
the temperature, the density of states, or even the system considered, increasing the number 



of quantum states. There are some interesting possible directions to 



present investigation, such as the study of non-Markovian systems 



Nakajima-Zwanzig projectors |26l l^ . 
tween the system and its environment 



ol 



ow to extend the 



19N24|. time-dependent 



ligher-order approximations for the interaction be- 



2a 



29j, and effects of other kinds of environmental 



noise. One problem which the present theory can not address is the possibility that the 
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measurement result be different from one of the eigenvalues of the meter |30| . 

The paper is structured as follows: In Sec. II we briefly present the proposed formalism 
and state the problem to be analyzed. In Sec. Ill, we present the approximate analytical 
solution to the master equation, followed by a derivation of its numerical counterpart in Sec. 
IV. In Sec. V, we discuss the results obtained with the expressions in Sec. Ill and Sec. IV 
and its physical implications. In Sec. VI, we present some perspectives for expanding the 
material that we have presented here and a conclusion. 

II. THE HYBRID MASTER EQUATION 

The Lindblad equation, 



lt PSB = ~h 



h, hs] + v (£fWf>* - \ {VpWf\ f lsB }) (1) 



where psB is the density operator describing the system and its environment, H is the total 

A ( s) 

Hamiltonian and the L - are the Lindblad operators, which act only on the system, is the 
most general form for a master equation fill \\\ . The second term on the right-hand side 
ofEq. (HJ, 



J 

is the Lindbladian operator acting on the density operator psb- The Liouvillian operator 



acting on p SB , 

i 



H, Psb , 

accounts for the unitary portion of the propagation, before the environmental degrees of 
freedom are traced out, and the Lindbladian represents the Markovian measurement dy- 
namics. 

We begin with a system S and its environment B, whose interaction is described by the 
Hamiltonian: 



Hsb = 2^ S k B k , (2) 
k 

where, for each index k, operates only on the system 5* and B k , only on the environment 
B. The form of the interaction, Eq. ([2]), is general enough, assumed by both amplitude 



and phase damping models 15|. Here, to account for non-Markovian noise, we choose to 
treat the environmental interaction as part of the total Hamiltonian, H, appearing in the 
Liouvillian term of Eq. ([!]). Thus, we write 



H = H B + H SB + H s , 

where H B is the environmental Hamiltonian and Hs is the system Hamiltonian. All the 
information read by the measuring apparatus, which is assumed to be Markovian, will be 
accounted for by the Lindbladian term of Eq. (CQ) and, as usual, the Lindblads Lj will 
act only on the Hilbert space of the system, since we are interested in measuring system 
observables only. Our aim is to obtain an equation for the time evolution of the reduced 
density matrix of the system, p$, that is, 



p s = Ti B {p SB } . 

For any density matrix operator X, let the superoperators B, S, and F be defined, respec- 
tively, as 



BX 



i 

h 



H B , X 



(3) 



SX 



i 

h 



f ^ ( Lf>XLf* --\L 3 



(S)tf(S) 



and 



FX 



x 

h 



Hsb,X 



It is easy to show that BS = SB and e st+Bt = e st e Bt = e Bt e st . 

Next, we will also use the Nakajima-Zwanzig projector superoperators 
defining action of the Nakajima-Zwanzig projector P is written as 
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(4) 



PX(t) = p B (t )®Tr B {x(t)} 



12| . The 



(5) 



for any X (t) and any initial time t - With these projectors and algebraic manipulations, we 
obtain the general equation, 

d - 



dt 



Pa (t) 



dt' 



PG (t) G (t') Pa (t) 



(6) 



where 

& (t) = e- k ~ k psB (t) (7) 

and 

6(t) = e -St-kpJt + k^ (g) 

Evidently, according to Eq. ([?]), once a (t) is calculated, ps (t) can be found by the action 
of e St on the reduced a (t) , that is, 

p s (t) = e^Ti B {a(t)}. 

As we explain in the Introduction, here we consider a two-state system interacting with 
a bath of harmonic oscillators, that emulates a phase- noisy environment. Thus, we take 

H s = hu cr z , 

H B = hJ2^kb{b k , (9) 



and phase- damping interaction 



15] . that is, 



Sk = ftcr z , 
B k = 9ki>l + g* k b k . 

where Uq and the uj k are real constants, b k and b k are the annihilation and creation bath 
operators, g k are complex coefficients. Here and next, we will use the Pauli matrices a z and 



1 \ | 1 



0-1 \ 1 



Hence, let us define the operator 



B = = E + aik) . (io) 



k k 

Therefore, the interaction can be written in the simplified form: 

Hsb = ha z B, 
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that is, 

H S b = h^a z (g k b{ + g* k b h 

k 

In the case of a finite temperature, we take the initial state of the environment as given by 

p B = j-Y[ e -'^i\ z B = nr^- (u) 

fc i 
III. SOLUTION OF THE HYBRID MASTER EQUATION 

In this section, we solve Eq. <Q analytically in the cases where l/ 5 ) = \a z at finite 
temperature, and = \a x , under the assumptions that H$ = and T = 0. The quantity 
Pa (t) appears on both sides of Eq. (JH]) and we can simplify it: 

Pa (t) = Pe~ §t - M p S B (t) = e~ §t p s (t) Tr B [e~ M p B j p B = e~ §t p s (t) p B . 
Now let us define the operator: 



R(t) = e- bt p s (t). (12) 



Hence, 



Pa (t) = R(t)p B . 

Therefore, to recover the reduced density operator of the system, we apply e st to R (t) : 

p s (t)=e St R(t). (13) 

An unusual aspect that should be clarified is the action of the superoperator exponentials, 
e st and e Bt . Let us consider, initially, the time-independent operators X' and X, related by 
the operation e Bt : 

X' = e k X. (14) 
When we take the time derivative of X' and use Eq. (|14|) . we obtain 
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Now, if we consider the definition of the superoperator B, Eq. ([3]), we obtain the elementary 
Liouville-von Neumann equation, that is, 

d 



X' — — — 
di h 



whose solution is easily determined for a time-independent Hb, as in the case of fl9]): 



X' = e Bt X = e-^Xe^. 



(15) 



A. Expanding the integrand that appears in the hybrid master equation 

In view of Eq. (fT2|) . the integrand can be written as PG (t) G (t') R (t) ps- From Eq. ([8]) 
we obtain 



G(t)G(t')R(t)p B = ie- §t e-^F\e § ^\(e §t 'R(t))a 



e m p B )B 



- ie~ §t e- k F {e § ^ [a z (e k 'R(t))]) 



x <^ e 



B ( e m 'p B 



(16) 
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From Eqs. (jlj) and (131) , we can rewrite Eq. (fT6l) as 



G (*)<?(*')£(*) P S 



e st 'R(t))a 



Tr B (e-^S le*^ 



x lr B 



Ps 



e 5t 'i? (t) ) *. 



cr 2 



x Tr B (e-« je^*"*') 



e^ps ) 5 



5 \ ®PB 



(II) 



e~ st a 7 e s ^ 



t z [e st 'R{t) 



x lr B 



Tr R je^-*') 



5 ( 



Pb 



Tim 



a z (e st 'R(t) 



cr T . 



x Tr R {e- Bt \e B ^ 



B ( e B % 



B \ ®p B . 



7m 



(17) 



It is interesting to note that, in Eq. ffTTj) . the actions of S and i? are completely sepa- 
rated in each term appearing on the right-hand side. This proves extremely valuable in the 
calculations that follow. 



B. Tracing out the environmental degrees of freedom 



For the sake of convenience, let us analyze, firstly, the environment quantities appearing 
in Eq. (|17|) . According to Appendix A, the partial trace over the environmental variables 
gives: 
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PG it) G it') Pa (t) 



e~ st aAe s ^ 



x 



e st 'R{t))a 
e k 'R(t))a~ 

^ \gi\ 2 jcoth (^p-^J cos [u)i (t - t')} + i sin [u x (t - t')} 



_ e -st e s(t-f) 



<Jz \ » Pb 



e -st I e S(t-t') 



,S'f 



(7, 



cr z e 



x 



^2 \9i\ 2 jcoth (^p-^J cos [wj (* - If)] -isin [u t (t - t')} 



C. Introducing a continuous density of states characterizing the environment 



In Eq. (1TB]) . if we adopt the general definition of the density of states as 
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then the sum over the index / can be replaced by an integral over a continuum of frequencies: 
PG it) G it') Pa it) = [ duJ (u) (coth f ^ J cos [u (t - t')) + % sin [u (t - t')) \ ® p B 



e St 'R(t))a z 



(A) 



e -st\ e s(t-t>) 



e st 'R (t) ) a 



(Jr. 



(B) 

+ I dujJ (cu) <{ coth 



h/3u 



cos [tu (t — t')} — i sin [u (t — t')] > ® pB 



x <e- st e^"*') 



a z [e st 'R(t) 



(C) 



a z [e st 'R(t) 




(20) 



(D) 

Here, to obtain an analytical solution, we choose the Ohmic density of states (fT9|) : 



J (oj) = rjcue' 



(21) 



where rj, u c ^ and 77 is the constant that gives the intensity of the coupling between the 
system and its environment. 

D. Reduced density operator describing the system 



To obtain the action of the operator e , it is necessary to solve Eq. ([I]) without the 
environment. Then, we take 

H = H s = ftu fr z 
and, in the Lindbladian, = Xa z or = \a r . 
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In the case of l/ 5 ) = Xa z , the solution is simple and can be found in Ref. 

P [ z ht) = P j?(o), 

p® (t) = pg (0) e-^e-^, 



(22) 



where the upper index (z) indicates that the solution is written in the eigenbasis of a z . 
For = Xa x , the solution is more complicated: 



+ 



p$ (0) cosh (yA 4 - Aoo^t 

(0) 7P=| sinh (V ^ - 

rP Sf (0) sinh (V A 4 - 



(23) 



To analyze the result of a measurement it is natural to represent the density operator in 
the eigenbasis of the measuring apparatus. Accordingly, in the present case, we use the 
eigenstates of a x , that is, 



_ !+)+!-) 
_ I+H-) 

The change of basis is performed with the eigenvectors matrix 




M 




M" 1 



(24) 
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and the result is: 

» 



i + e -*t j C o S h (^A 4 - 4u$t) Re {pg (*)} 
sinh (v^4^ft) Im {pg (t)} 



- v /A 4 -4w ( 



= sinh (V A 4 - 4w 2 t) Re jpg (i) 



(25) 



Pi? (*) 



2 Pil ) W~ 1 c -2A 2 t _ - e -A 2 



* {cosh (yA 4 -4w 2 t) Im jpg (t) 
v/ ^smh( v ^4^t)Re{pi?(t)} 



2up 



= sinh(VA 4 -4o; 2 t) Im{pg(t)U. 



V A4 - 4 ^ 



E. The case of = Xa z 

Here, we consider the evolution of the system in contact with a thermal reservoir at 
arbitrary temperature, i.e., we assume that the initial condition of the environment is given 
byEq. JHJ. 

Let us write 



R(t) 



(26) 



Ru R\2 
R21 R22 

where, for notational convenience, we take Rij = Rij (t) . Then, using Eq. (l22l) . we obtain 



e st 'R(t) 



R11 R 12 e 2xH ' e-^ 1 ' 

?2ie 2AV e ^ t' ^ 

Substituting this into Eq. fl20l) and manipulating the terms according to Appendix B gives 



d ( R11 R12 
dt 1 R21 R22 




J* dt' duJ (to) cos [u it - t')] coth (^) ■ (27) 



1. Obtaining the populations 



According to Eq. (127|) . the populations are independent of J (u) and can be immediately 
evaluated, giving 

13 



d 
Tit 



Ru = R u (t) = Ru (0) 



where % = 1,2. Then, using the constraint that the trace of the density operator must equal 
unity, we obtain 



Pn(t) =Pn(0), 
p 22 (t) =1- Pll (0). 



(28) 



2. Obtaining the coherences 



From Eq. (|2yI) it follows that the non-diagonal elements satisfy 



d_ 
dt 



Rij = —Ar/Rij (t) / dt' / ducue~~ cos [u (t — t')} coth I — 
Jo Jo V 2 

where i, j = 1,2 and % ^ j. From the procedure of Appendix C, we obtain 



(29) 



Pu{t) =Pl2(0) 



r (^+^) r (^-^) r (^+ 1 +^) r (^+ 1 -^) 



2»7 



e -2X 2 t e i2ui t 



(30) 



where T denotes the gamma function. Simplifying Eq. (|30|) . according to Appendix D, gives 



P12 (t) = P12 (0) 



(u c ty 



n=l 



1 + nu c f3h 
+ n<jj c (3h + cj c t; 



2JJ 



(31) 



The consistency of Eqs. ( !30|) and ( )3~Tj) with the case without environment can be verified 
by noticing that, in the limit rj — > (i.e., in the absence of environment), they reduce to Eq. 
as expected. 



F. The case of l/ 5 ) = Xa 7 



From Eq. (12"3"1) the action of e st on an operator 



x 



x u x 12 



%21 %22 



gives 



14 



where 



e St x 



Sll S12 



■521 ^22 



11 2 2' 

_ 3^11-^22 ^~2A 2 t i XU+X22 

$22 — n — e H 5 — , 



(32) 



with 



Sl2 
S21 



[fi cosh (fit) X12 — i2a;o sinh (fit) x 12 + A 2 sinh (fit) a? 2 i] 
[fi cosh (fit) X21 + i2uo sinh (fit) x 2 x + A 2 sinh (fit) xx 2 ] 



(33) 



fi = a/A 4 -Aul 



According to the procedure explained in Appendix E, using Eqs. (152"]) and (133]) in the 
case of the operator 



R{t) 



results in the final differential equation: 



ru (t) r 12 (t) 
r 2 i (t) r 22 (t) 



(34) 



dt 



rn(t) r 12 (t) 
r 2 i (t) r 22 (t) 



jT dt' jf rfuJ (w) cos [w (t - t')] coth f^r^l 



x 







Q* 2 {t,t')r 12 {t) + Qt {t,t')r 21 (t) 



Qi (t,t')r 12 (t) + Q 2 {t,t')r 21 (t) 







(3 



where 



Qx (t, t') = A a it) [K{ (t - t') A7 (f) - A 2 (t - t') A 2 (t')] 
+A 2 (t) [A 2 (t - t') A* (*') - Ax (t - if) K 2 {t')\ 

Q 2 (t, t') = K x (t) [A* (t - t') A 2 (f) - A 2 (t - t') Kx {if)] 
+A 2 (t) [A 2 (t - t') A 2 (f) - Kx (t - tO Kx {if)] 



(36) 
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and 

K\ (t) = Q cosh (fit) + i2u sinh (fit) , 
K 2 (t) = A 2 sinh (fit) . 

1. The populations represented in the eigenbasis of a z 

Analogously to the case of L^ s ' = \cr z and T ^ 0, the populations do not depend on 
J (uj). From Eq. fl35l) we obtain 

^-ru = ^ va (t) = r u (0) , 
where i = 1,2. We notice that ru (0) = pf} (0) and, using Eq. (fT3|) . we can write: 

f& {f) = 2^-1^, + 1 

5. T/ie coherences represented in the eigenbasis of a z , for T = and ujq = 

Since ri2 (i) is the complex conjugate of r 2 \ (£), we only need to calculate one of them. 
Let us introduce the new variable r = t — t'. Hence, it follows from Eq. ( 135]) that 



d A T) f , f°° , _«i , , , /0fe«A 

77^12 (*) = -4ttt / dr \ duue ^ cos wr coth — — 

M ^ JO JO V 2 / 

x [Q 1 (t,t-r)r 12 (t) + Q 2 (t,t-r)r 21 (t)\. (37) 



There is no analytic solution for this equation at a finite temperature. However, as detailed 
in Appendix F, we have found the following result for T = and uq = : 

+ zlm U$ ^ e ^ e ^ 90 t e ~A^[A + (t) + B +{ t)]^ (3g) 
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where 

A + (t) 



„2A 2 t 



'gi (f) df 



+ 



l 

4A 2 

1 
2A 2 



Re i e 2 u c 



Re { e 12 



e 4A2i r I 0, 2A*t + i2— I - e~ l ^r I 0, -2A z t - z2— 



A 5 



*4^ n 



A 2 / A 2 
i2— r [ 0, -2A 2 t -i2— 



A- 



UJ r 



r 1, -2A 2 t - i2- 



A- 



+ tRe <( e" i2 ^r I 0, -2A z t - i2— I }> + Re <J e~"^ C i + e lz ~c 3 



uj, 



-i2±- 



\2 

J2— 



A- (t) 



-2\ 2 t 



9i (?) dt' 



4A 2 



Re { e i2 "c 



-4A 2 i 



2A 2 



Re<U 



i2- 



A 2 



A 2 

r ( 0,-2A^ -i2— 
A 2 



e i4 £r I 0,2A z t + z2— 



A 2 
>_ 



%2—T 0, 2\H + i2— - T 1, 2A^t + %2— 



UJ, 



A- 



ijj, 



+ tRe<J e i2 ^r I 0,2A^ + z2 



A 5 



UJ, 



i i2 — * . -i2 — 

Re ^ e *Cj + e *c 3 /, 



5 + (t) 



„2A 2 t 



'^2 (0 eft' 



o 

1 

4A 2 
1 

2A 2 



Re^ e 



i2- 



A M-l,2A 2 t + 22-^1 +e~ i4 ^r 



Co', 



-l,-2A 2 t-z2- 



A- 



Re<e 



-i2- 



\2 / \2\ / \2 

i2— r -1, -2A 2 t - %2— + T 0, -2A 2 t - 22— 



0J r 



ijj, 



UJ, 



- tKele 



-i2- 



-l,-2A 2 t-z2- 



uj, 



T-) i — «2— iz — 

— Re < e "c C2 — e * C4 



/2- 



= f e- 2XH 'g 2 {t')dt' 
Jo 



'0 
1 

4A 2 



Re -4 e j2 ^c 



t2- 



A 2 
1 



e -4A»t r _ 1; _ 2A 2 t _ z2 — + e l4 ^r -1, 2A 2 t + i2 



A 2 



A- 



A 5 



i2— r -1, 2A^ + z2— - T 0, 2XH + i2— 



A 5 



tRe { e i2 ^T I -l,2A z t + i2- 



A^ 



1 i2^- * -i2^ 

Re < e <"= Cn — e u = c, 



Co' 



4 r ' 



and 
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Cl 
C4 



1 

'2A 2 
1 

~2A 2 
1 

'4A 2 
1 

'4A 2 



2^r(o. 

wo 



U>C 



—I 



2A 2 
2A 2 



r i 



;2A^ 

(0,-z" 



r o,i^ 

UJ C 



^r(o,-i^ 



r 



;2A 2 



2A 2 



3. The density matrix elements represented in the eigenbasis of a x , for T = and = 

Let us change the previous result for the coherences to the eigenbasis of a x . Using Eq. 
(125]) with wo = we obtain 

P<$ (t)= \ + Re {pSJ (0)) e -^>?9ot e ^[A-{t)-B-{t)]^ 



2pg(0)-l 2A 2 f 

O fc- 



-zIm<|pSJ(0)^ e 



-2A 2 t p 87?A 2 g tp-4 ? ?A 2 [ J 4 + (t)+B + (t)] 



(39) 



IV. NUMERICAL METHOD 



The approximate analytical results detailed above have been compared against the exact 
numerical solution of the Lindblad equation, Eq. ([T|). The algorithm we use derives from the 



25| . which is adequate for numerically solving a first-order 



superoperator-splitting method 
differential equation of the form: 

^-PSB (t) = Ap SB (t) + bp SB (t) . 

at 

As long as the two superoperators A and B are time-independent, the solution for the 
differential equation is 

Psb (t) = e( A+B yp S B (0) . 

In case this exponential superoperator cannot be analytically found, but the e At and e Bt 
can, we may use the approximation 



A+B)At 



e AAt £ BAt 



O (At 2 ) 



to expand the solution in terms of a product of N short time steps of length At. The alternate 
application of these two superoperators, e AAt and e BAt , comprises the superoperator-splitting 
method. 
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For a finite-time measurement, the superoperator A represents the Liouvillian and B 
represents the Lindbladian. The form of the exponential of the time-independent Liouvillian 
is straightforward, while e BAt behaves in the way described by Eqs. fl32|) and fl33|) . As the 
behavior of both populations in the a z eigenbasis can be exactly described by a simple 
analytical formula, only the coherences need to be calculated numerically. Therefore, all 
the relevant information about the density matrix can be represented by a pseudo-spinor 
containing both coherences, on which the alternate action of e AAt and e BAt will be equivalent 
to matrix products: 



e 



K+i (At) A +1 (At) + K- X (At) A_ a (At) , (40) 



P21 (NAt) J 1 J \ p 2 i (0) 

where 

K (At) X = e ~ i ^k^k(h+qgk/u k )\b k +qg k /u k )Atj£ e iY, k Uk(h-qgk/u k )\h-<]gk/u k )At 
and we defined the square matrices A± (At) as 



A +1 (At) 
A_i (At) 



with 



b q (At) = \ {f M + qe 



6+i (At) 6_i (At) 




6_i(At) 6+i (At) 

X 2 At , „-X 2 At 



2 

The binomial in Eq. (|40j) may be expanded in two-term summations. Also, as we are 
not interested in the total density matrix, but in the reduced part referring to the two-state 
system, we may take the partial trace over the degrees of freedom of the environment, to 
find: 



P12 (NAt) 
P21 (NAt) 



N 



X 



g— A 2 (NAt) £ ... £ Y[[A qn (At)] 

gie{-l,l} q N e{-l,l}n=l 

^IswHt!)' (41) 
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where we are considering the system initially separable from the environment, which starts 
in the vacuum state (T = 0). 

Employing the result for the trace given in Appendix G and the easily-verified matrix 
identity, 

A r (At) A s (At) = b rs (At)(a x f~ rs)/2 A s (At), 

we find the final form for the algorithm, which is comprised of a sum of 2 N terms, each 
containing N 2 (N — 1) factors and a product of up to two matrices: 



P12 (NAt) 
P21 (NAt) 



e -X 2 (NAt) ... ^ 

^{-1,1} 

N ( of a^- 2 , d o i i 2 A ~ m 

t-j- I ^ _|_ ^ {cj c At) + (1 — 2 \m — n\ ) 



21 2 



[(w c At) + |m — to | ] 

n m ^f- qiqN)/2 a n (At) [ pi2 f \ , (42) 

V P21 (o) 



m,n=l 
N-l 

X 

71=1 



which can be implemented in any programming language capable of handling large floating- 
point numbers and of calculating powers and exponentials. The results are then translated 
into elements of the density matrix in the basis of eigenvectors of a x by a simple transfor- 
mation, Eq. (124)) . 



V. RESULTS AND DISCUSSION 



The final analytical expression for the finite-time a x measurement given in Eq. (1391) is 
expected to agree with the numerical solution Eq. (1421) for weak coupling with the environ- 
ment, in our case represented by the dimensionless constant r]. This fact is confirmed in 
Fig. 1, which compares the time evolution of the population p\\ it) using both methods for 
different values of the system-environment coupling constant. There it can be seen that the 
discrepancy between the two graphs grows with rj. 

Having confirmed the reliability of our results, we proceed to compare the noisy measure- 
ment described by our formalism with the case where the system is not being measured. As 
can be seen in Fig. 2, systems evolving only under the influence of the environment (A 2 = 0) 
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0.86 




0.2 0.4 0.6 0.8 1 

t (units of a> c " 1 ) 

Figure 1: (Color online) Comparison between the time evolution of the population in the eigenbasis 
of a x , (t) , given by the analytical (crosses) and numerical (bullets) methods for A 2 = 4a; c and 
different values of r/. The system is initially in the pure state (|1) + e J7r//4 |2)). 

always have faster rates of change of their populations than those which are also interacting 
with the measurement apparatus (A 2 > 0). Moreover, the stronger the measurement, the 
slower the rate of change of the population. Therefore, a measurement of the observable a x 
helps to maintain the population closer to its original value. 

The meaning of this phenomenon is straightforward: a difference between the initial 
value of the population (instant to) and its value at the end of the measurement (instant 
tf) represents the probability that the measurement will give a wrong result. This is so 
because, if the initial value pn (to) and the final result pn (tf) are different, then there will 
be a fraction of e = \pu (tf) — pu (to) \ systems in the statistical ensemble that present final 
collapsed states different from those that would be obtained if we had measured them with 
an ideal instantaneous process at instant to- We shall henceforth refer to the difference e as 
the error. 

The fact that a finite-time measurement helps to preserve the initial value of the popula- 
tion of a system also shows that a naive approach to a noisy measurement will overestimate 
the error. Given that a measurement ends at a time tf and that the system has been inter- 
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Analytical results 
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0.4 0.6 
t (units of eo c ) 



Figure 2: (Color online) Time evolution of the population in the eigenbasis of a x when there is 
an external environment causing phase noise (crosses for the analytical method and bullets for 
the numerical). The dashed line represents the evolution of a system that is only subject to the 
interaction with the environment (A = 0). In all graphs, the evolution of the system that is not 
subject to measurement diverges faster from the initial value than those systems that are being 
measured. The initial conditions are the same applied to Fig. 1. 
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acting with the environment since instant to < tf, modeling the noisy measurement process 
as a period of interaction with the environment followed by an instantaneous measurement 
would lead to an error much larger than that obtained by considering that the system is 
measured continuously during a period tf — to- Therefore, a description of the measurement 
as a continuous process is essential to have good estimates of the error. 

Finally, it is important to notice that this protection against error depends on whether 
the observable measured commutes with the interaction Hamiltonian. In the measurement 
described above, the observable a x anti-commutes with the operator a z in the Hamiltonian. 
But when we measure an observable that does commute, the populations remain constant, 
and the only effect of the measuring apparatus is to strengthen the effects of the decoherence 
caused by the environment, as can be seen from Fig. 3. 

In short, when the observable measured is a z , the measurement apparatus reinforces the 
effect of the environment, thus leading to faster decoherence. When the observable is a x , 
the two effects compete against each other, leading to a smaller error than that expected for 
a system that is not subject to measurement. 

VI. CONCLUSIONS AND PERSPECTIVES 

In this work, we have analyzed a two-state system subject to finite-time measurements 
that commute or anti-commute with the interaction Hamiltonian which couples the system 
with a phase-noisy environment. In the first case, we have shown that the complete analytical 
results for a finite-time measurement in any temperature - seen in Eqs. ff28|) , (130)) and fl3Tj) - 
give only faster rates of decoherence as the strength of the measurement grows. In the latter 
case, we have arrived at an analytical expression for zero temperature in Eq. fl39l) which, 
together with the numerical method in Eq. (H2l) . has made it possible to conclude that 
the initial state of the system under consideration is protected by the measurement. This 
result, a consequence of the fact that two simultaneous measurements of non-commuting 
observables interfere with each other, allows a deeper understanding of the measurement 
process by opening some interesting perspectives in the area of protection against errors and 
demonstrating that it is necessary to take into account that the measurement is a continuous 
process rather than an instantaneous collapse in order to have reliable estimates of the error. 

Our approach of considering the system interacting with the environment, but the ap- 
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Figure 3: (Color online) Time evolution of the absolute value of the coherence in the eigenbasis of 
u z of a system suffering phase noise from an Ohmic environment (77 = 0.05), while the observable 
a z is continuously measured. The curves were calculated using the analytical result of Eq. (|31D . 
The more the strength of measurement increases, the faster the coherence goes to zero. The initial 
conditions are the same applied in Fig. 1. 
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paratus interacting only with the system, through our treatment of the Lindblad equation, 
is novel. However, we have treated the problem with a series of simplifications, namely: 
zero temperature for the initial state of the environment, Ohmic density of states, and the 
neglect of the system Hamiltonian when measuring an observable that does not commute 
with it. Subsequent studies may use our method in more general cases, even varying the 
number of states of the system considered. A first step in the direction of studying sys- 
tems interacting with environments initially at temperature T > was given in Eq. fl37|) . 
which shows the differential equations for the case of finite temperature that may serve for 
numerical resolutions that could reveal new effects. 

The present formalism does not allow us to approach the possibility of the state mea- 
sured being different from an eigenvalue of the system, as considered in some works on 
quantum theory of measurement {3^. However, it is still open the possibility of expand- 
ing our method for non-Markovian systems, time-dependent Nakajima-Zwanzig operators, 



and to analyze the quantum Zeno effect 



33 



-35|. In particular, the quantum Zeno effect, 



viewed as a dynamical phenomenon |36|, is closely connected with the fact that a finite-time 
measurement protects the initial state in the case of a measured observable that does not 
commute with the perturbing Hamiltonian. Roughly speaking, a finite-time measurement 
could be thought of as an approximation for a sequence of repeated instantaneous measure- 



ments, as in the bang-bang decoupling method [371 ]. which has been generalized to arbitrarily 
fast and strong pulse sequences, requiring no symmetry. Therefore, we expect to be able 
to connect the findings of Ref. js7| with our present approach, clarifying how a general, 
finite-time measurement dynamically induces the quantum Zeno effect to protect the initial 
state being measured. 



Acknowledgments 

C. A. Brasil acknowledges support from Coordenagao de Aperfeigoamento de Pessoal de 
Ni'vel Superior (CAPES), Brazil. 

L. A. de Castro acknowledges support from Fundagao de Amparo a Pesquisa do Estado 
de Sao Paulo (FAPESP), Brazil, project number 2009/12460-0. 

R. d. J. Napolitano acknowledges support from Conselho Nacional de Desenvolvimento 
Cientffico e Tecnologico (CNPq), Brazil. 



25 



Appendix A: Calculation of the environment degrees of freedom (Sec. iii-b) 
Using Eq. (fT5|) . we can write: 



(/) Tr B \e-^B \e 6 ^ 



e Bt 'p B ) B 



Tr B { e^Be-^pse^Be-^' (43) 



(//) Tr B {e~ k {e i (*-*') [(e k ' p B ) B 



B 



Tr B <^ e J r '.Be fi 7? s e l " .Be 



ft 



44) 



(///) Tr B {e^B j e V*') (e^ Pb)] }} 



Tr B <! e 1 ^' Be-^' pBe^Be-^ 



45) 



(IT/) Tr B (e-^ {e^" 4 ') 



5 e 



Tr B {e^Be^/be'V'&f 



>46) 



Let us insert, then, Eqs. 
terms: 



(l45p . and (146]) into Eq. (fT7|) and group the similar 



PG (t) G (f) Pa (t) 



e st a 7 \ e 3 ^ 



e St 'R(t))a 



_ e -st I e s(t-t>) 



e st 'R (t) ) a 



0-, 



x Tr B \ e^Be-^pBe^Be-^*' } ® p B 



St J p S(t-f ) 



(7- 



o, e 



^ (*))]}} 



x Tr B \ e l ^ v Be~ l ^ v pBe'^Be- 1 — 1 



Pb- 



(47) 
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The two terms in Eq. (1471) involving the trace over the environmental degrees of freedom 
have the same form. Then, it is sufficient to evaluate the first of them and perform the 
change 

t -»■ a 

a -> t 

to obtain the second. To take the partial trace, we use the Fock-state basis \m) <8> ^2) ••• 
]1 \n m ). Thus, 



and we obtain 



3l| 



Ik 



\ \n m ) 



Y,(g$ k + 9kbl)e < 1 = J2 [Ae-^ + gkble**) . (48) 

k k 

if we specify Hb and i? as given by Eqs. (Q and (fTUj) . By Eq. (JUJ), the trace over the 

environmental degrees of freedom becomes 



Tt b < e ft *5e « p B e fi -Be ft 



E 



ni,n 2 , 



x 



J |n m ) 



With (TTll. we will obtain 



Tr B < e ft .Be » p B e ft .Be 1 » 



E 
E 



.9/1 



.9/1 



e i0Jl(t-t') _|_ ( 



iui{t-t') 



-iu>i(t-t')- 



e iudt-t) + 



2cos[c^ (t -t')]" 

e /i/3w, _ J 



or, more simply, 



Tr B <| e^e^he^Be^ \ = £ l^l' {coth cos [a* (t - t')] 

+ zsin[wi(t-0]}- (49) 
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As explained above, the change of varibles 

t ->■ a 
t' -»■ t 

in Eq. (|49|) gives us the other necessary term for Eq. (T4"Tj) . that is, 



jilE.*/ A -i2M.fi ~ ;ia f A -i2Mf 

Tr B ^ e h fie ft p B e « fie » 



= X] I 5 ''' 2 { coth (~Y^) cos ^ ~~ ^ 
- zsin[^(t-t')]}. (50) 

The final result, Eq. ([18]) . is achieved by inserting Eqs. (|4~9]) and (150]) into Eq. (|47|) . 
Appendix B: the case of L( s ) = A<7 2 (Sec. III-E) 

Using Eq. (125]) we can evaluate the terms (A), (fi), (C), and (fi) of Eq. (120]) . Therefore, 



(A) e-% 2 e^) (e^'fi(t))<j 



fill — fil2 
— fi21 fi22 



(51) 



(fi) e- 



St J pS(t-t') 



e Si 'fi (t) ) a 



(7 - 



fill fi 
i?21 fi22 



12 



(52) 



(C) e 



-st J e S(t-t') 



cr 2 



,St' 



fi(t) 



o. 




(53) 



and 



(fi) e 



a z e 



-,st 



'fi(t) 



fill fil2 
fi2l i?22 



(54) 
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With Eqs. ([SI]), d52J, (JS3J, and ([54}, we can simplify Eq. ([20]): 



PG (t) G (f) Pa (t) = -4 I duo J (co) | ^° ) cos [w (t - t')} coth ) <g> p^55) 



Hence, we substitute Eq. f[55l) into the original Eq. ([6]) and obtain Eq. ([2 



Appendix C: the coherence for the case of = Xa z (Sec. III-E.2) 
In ( !29|) . for the sake of convenience, we perform the change of variable 

r = t-t', 
dr = —dt' 



with 



that gives 



dt' 



rt 

dr = dr, 
Jo 



d C f°° 

-j-Rij = —4:T)Rij (t) ^ T dcocoe~^ cos (cot) coth 
dt J J 

The first step is to solve the time integral, that leads us to 



/3hio 



d n „ , . f°° , i \ , fPhu 

-j-Rij = —QrjRij (t) / due w = sin (cot) coth I 

at jo \ 2 



The frequency integral, evaluated with the help of Ref. 



32j, gives 



where 



d 



i2-^Rij (t) 
1 



dt' 11 ' ph' 



co c (3h f3h 



^ ) 



t 



+ V 



co c (3h (3h 



t 



ijj(z) = ^-\n[T(z)} ,Re{^} > 0. 
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The solution is, therefore, 



Ru(t)= R 12 {0) 

R2l(t)= i?2l(0) 



Xi^+^) r (^-^) r (^+ 1 +^) r (^+ 1 -^) 



r2 (^ + l) 



r2 ( w c /3fi> 



2,; 



1 2J7 



(56) 



The coherences are obtained by performing the transformation of Eq. (1131) to Eq. pb]). 
Since p 2 i (i) — Pi 2 (0> we consider only the element p 12 (i) and the result is Eq. (130]) . 



Appendix D: Simplifying the expression of the coherence for the case of = A(7 2 
(Sec. III-E.2) 



Equation (130]) can be simplified by using the properties of the gamma function 32]. With 
r (x + 1) = xT (x), we obtain 



Pl2 (*) = Pl2 (0) 



2/1 , • t \ T 2 I 1 „■ t 



oj c/ 9^ ' phi \ oj c /3h 0h 



F 4 f _JL_ 

1 u> c f3H 



2') 



2\ 2 tJ2u t 



Now, using the product representation 



rW=e- c, in& Re{4>0, 
fc=i fc 



where C is the Euler constant, 



C= lim 



E - - In (s) 

m=l 



0.577215..., 



we see that 



r2 (^T8ft) 



n=l 



resulting, finally, in Eq. MM, which is suitable for numerical implementations. 



Appendix E: the differential equations for the case of l/ 5 ) = Xa x (Sec. III-F) 

Using the definitions (132]) and (133]) over (1M)) . we have 



30 



where 



and 



,(o) 

12 



e st 'R(t) = e st ' 









1 rg> 




\ r 2 i 


^22 j 


H 







».( ) _ ru-r 2 2 „-2\ 2 t' I rn+r 2 2 

'XI — 2 C 2 ' 

~(°) _ r -2X 2 t' _i_ r-n+r 2 2 

22 — 2 C "r" 2 ' 



[fi cosh (fit') r 12 - i2w sinh (fit') r i2 + A 2 sinh (fit') r 2 i] 



r g) = £-11 [fi cosh (fit') r 2 i + z2w sinh (fit') r 2 i + A 2 sinh (fit') r 12 ] 
Now we can evaluate the terms between brackets on fl20l) . that will furnish 



(57) 



(58) 



a z [e st 'R(t) 



'21 '22 



and 



J(t-t>) 



e st 'R(t))a 



'21 '22 



where, using (|57j) and (158)) . 



n 

22 



(0) , (0) n , (0) 



2 

,(°) 



-e 

,(°> 



+ 



'n + r 22 ; c -2X 2 (t-t') _|_ rff- 



,.(<>) 

'22 



2i _ n 



(2) _ r^+rff 2A2 q- n rj°^ 
'11 _ 2 C ^ 2 

(o),^(o) „,„„ , (N Jo) (o) 



,(2) _ r^+r 22 



22 — 2 e ' T 2 



2\ 2 {t-t') , ~ r 22 
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and 



„(2) 
12 

„(2) 
21 



-A 2 (t-i') 



n 

-\ 2 (t-t') 
n 



{fi cosh [fi (t - t')] r { $ - i2u sinh [fi (t - t')] rj? - A 2 sinh [ft (t - t')] rf? } 
| ft cosh [ft (t - if)] r|? + i2u sinh [ft (t - If)] r|? - A 2 sinh [ft (t - if)] rg } | , 



i.e., 



n 

22 



„(2) 
11 

„(2) 
22 



and 



12 
21 



„(2) 
12 

„(2)' 
21 



The progressive calculus of (A), (B), (C) and (£)), then, will furnish: 



(^4) 



.st 



R(t))a 



/y-i /in (^) 

'21 '22 



(59) 



with 



,(3) 
11 

r (3) 
'22 



r ( 2 ) +r < 2 ) „, 2 , r ( 2 )_ r ( 2 ) 
r ll + r 22 p 2A 2 f _|_ r ll r 22 



r ( 2 ) +r ,( 2 ) ox2 + 
r ll + r 22 ^2A 2 t 



11 T 22 



and 



.(3) 
12 

,(3) 
21 



(2) 
21 



ft cosh (fit) r{J + i2w sinh (fit) rjj + A 2 sinh (fit) r, 
ft cosh (fit) r|? - i2w sinh (fit) + A 2 sinh (fit) r 



(2) 
12 



(B) 



e -st I e S(t-f ) 



R (t) a. 



(T- 



r 



r 



(4) _(4) 

11 ' 12 

(4) _(4) 

21 '22 



(60) 



with 



n 

„(4) 
22 



(2) 
11 



r ( 2 ) 
'22 



2 

„(2) 



+r£> 2 t + 



,( 2 ) 
li 



r ,2) 
' 22 



and 



,(4) 
12 



r (4) 
'21 



fi cosh (fit) r j 2 2 } 



i2o;o sinh (fit) r ^ + A 2 sinh (fit) r 21 
fi cosh (fit) r$ - i2oo sinh (fit) rg } + A 2 sinh (fit) rj? 



5 A ' 

n 
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We can observe that 



(C) 



with 



and 



(3) 
11 


= rg> 

and < 




_ r (4) 
— ' 12 


(3) 
22 


_ r (4) 
_ '22 


4? 


_ r (4) 
_ '21 



„(5) 
11 

.(5) 
22 



r (i) 

'22 



+ 



2 



'21 '22 



r (1) 
'22 



(61) 



,(5) 
12 

„(5) 
21 



n 



fi cosh (fit) rg + i2u sinh (fit) rJJ + A 2 sinh (fit) 
fi cosh (fit) r$jf - i2uj sinh (fit) r£ } + A 2 sinh (fit) rg 



with 



and 



(D) 



e - st a z I e s ^ 



„(6) 
11 

J6) 
22 



,St' 



R(t) 



"n+'-22 c 2A 2 t 



2 



'21 '22 



!2_g2A 2 t _|_ 



(62) 



,(6) 
12 

>) 
21 



IT 



(1) 

21 



"TT 



fi cosh (fit) r$ + i2w sinh (fit) + A 2 sinh (fit) r 
fi cosh (fit) r$ - i2w sinh (fit) r$ + A 2 sinh (fit) r 



(5) 
11 


= rj? 

and < 


>> 


_ r (6) 
— ' 12 


(5) 
22 


_ r (6) 

_ '22 




_ r (6) 
_ '21 



Again, we have another similarity, 

By using fl59|) . fl60|) . f lSTj) and f )62|) . with the similarity relations, we have then the equation 

dSSD. 
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Appendix F: the coherences in the case of = Xa x represented in the eigenbasis 
of a z , for T = and uj = (Sec. III-F.2) 

For T = and Uq = 0, the equation for the coherences becomes 
d 



with 



and 



dt 



r 12 (t) = A-q [sinh (X 2 t) I 2 (t) + cosh (X 2 t) I 3 (t)] r 21 (t) 
- Arj [cosh (X 2 t) I 2 (t) + sinh (X 2 t) I 3 (t)] r 12 (t) 



/ 2 (t) = 2A 2 sinh (A 2 t) # + A 2 sinh (X 2 t) g x (t) + A 2 cosh (A 2 t) g 2 (t) , 
J 3 (t) = -2A 2 cosh (A 2 t) # + A 2 cosh (A 2 t) g x (t) + A 2 sinh (A 2 t) g 2 (t) 



(63) 



(64) 
(65) 



, 2A 2 \/ 2A 2 
g = Re <^ exp z T O, z 



and 



2A 2 \ / 2A 2 

g 1 (t) = Re { exp ( 2A 2 t + % T O, 2A 2 t + % 

u c J V u c 



+ exp 



2A 2 t + z 



2A^ 



0, - 2A 2 t + z- 



,2A 2 



92 (t) = Re 



exp 









|exp 




r 









1, 2A 2 t + z- 



2A 



2\ 1 



2X 2 t + i- 



.2A 2 

u c 



u c 

>2 2A 2 
■1,- 2X 2 t + i 

u c 



Let us write the complex conjugate of Eq. (163]) : 



-r 21 (t) = 477 [sinh (A 2 t) J 2 (t) + cosh (A 2 t) J 3 (t)] r 12 (t) 

- Arj [cosh (A 2 t) J 2 (t) + sinh (A 2 t) I 3 (t)] r 21 (t) . (66) 

Adding Eqs. fl63l) and f l66l) and subtracting Eq. f l66l) from fl63l) . we obtain the decoupled 
system of equations: 
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|Re {r 12 (t)} = A V [sinh (\H) - cosh (A 2 t)] [J 2 (*) - J 3 (t)] Re {r 12 (t)} , 
|lm {r 12 (t)} = -477 [sinh (X 2 t) + cosh (A 2 t)] [J 2 (t) + J 3 (t)] Im {r 12 (t)} 
Simplifying Eqs. (164)) and (165]) gives 

|Re {r 12 (*)} = -4r/A 2 [2^ - e" 2 ^ (t) + e~ 2A2 ^ 2 (t)] Re {r 12 (i)} , 
|lm {r 12 (t)} = 4r/A 2 [2^0 - e 2xH gx (t) - e 2x2t g 2 (t)] Im {r 12 (i)} , 
whose solutions, in terms of the initial instant to — 0, are given by 



'Re {r 12 (t)} = Re jpg (0)} e' 8 ^^ eX p {47/A 2 jj e~ 2xH ' g x it') dt'} 
x exp |-4r/A 2 /„* e" 2 ^*'^ (*') di 



Im {r 12 (t)} = Im jpg (0)} e 8 ^' exp |-4r/A 2 jj e 2 ^! (f ) dt'} 
:exp|-4r/A 2 / *e 2A2 *'(7 2 (0 ^} ■ 



Calculating the integrals 



32j results in Eq. 



Appendix G: Calculation of the trace in Eq. ()41|) 

A convenient method for calculating the partial trace in Eq. (|4ip involves coherent states. 
For any sequence of q\, . . . , q^, this partial trace takes the form of a matrix element of the 
vacuum: 



N 



N-1 



.71=1 



T±b{] \K qn m |0)(0| = ; J <0| fc J e 



k n=0 
N 



,iu) k (b k —q N - n g k /(v k ) (b k -q N - n g k /ui k )At 



n=l 



-iw k {b k +q n g k /u k ) (b k +q n g k /uj k )At 



10), • 



The exponential operators present in this matrix element are best represented by displace- 
ment operators (gk/oJk) from quantum optics, which allow us to rewrite the partial trace 



as: 
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TV 



TV-1 



Tr B \] \k qn (At) |0)(0| = ] I <0| fc ] 



n=l 



k 

TV 

n 

n=l L 



n=0 L 



t L 9k 



9k 



9k_ 



SIS,, At ; 



-icu k b' b k 



|o>*. 



Any sequence of these operators applied to a coherent state \a k ) yields: 



D t ( gjh \ e -W&A*£ / 9h 



where we have discarded the complex phase factors due to the displacement operators. 
Repeating the procedure N times, we find that the partial trace is simply the inner product 
of coherent states: 



TV 



k \n=l 



n=l 



-inui^At 



(1 



Ju) k At 



\ 9k \ 
) QN-n+1 / • 



Such an inner product results in the exponential: 



N 



TV TV 



n=l 



Tib I ] K qn (At) |0> (0| \ = exp \ -8 ^ ]T g m g, 



JM 



sin 



m=l n=l 
X COS [(?Ti — 71) UlAt\\ 



where we have taken the limit to a continuous spectrum of frequencies, applying the defini- 
tion of the spectral-density function given in Eq. (TT9I) . 

For Ohmic spectral densities, Eq. (}2Tj) . the partial trace becomes the following exponen- 
tial of an integral: 



TV 



.n=l 



Tr B ] \k qn (At) |0> 



TV TV 

= J exp < 

m=l n=l 
X cos [(pi — Tl) id At] } . 



duj- 



sin 



uAt 



The remaining integral in dco may be solved by first noticing that f Q dr sin 
^ sin 2 This leads to a double integral that is easily solved as: 



LOT 



3G 



rAt poo 

4 / dr duo e _aj/Wc sin (ur) cos [(m - n) uAt] = In < H 

Jo Jo { [l 



2 (tu c Aty 



+ 



+ (m - n) 2 (w c At) 2 ] 
[l-2(m-n) 2 ] (w c At) 4 l 
[1 + (m - nf (u c At) 2 ] 2 J ' 



which gives the following result for the trace: 



I tt r s , A J , w , tt tt I 2 (u c Aty 2 + [1 - 2 (m - nf] 



21 ^ —VqmQn 



n=l 



m=l n=l 
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